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ABSTRACT 

In spectra of many Seyfert galaxies there is a wide emission line of Fe Ka- The line 
profile with two maxima supposes that the line emerges in innermost regions of an 
accretion disk around a black hole, hence, it is necessary to take into account General 
Relativity (GR) effects. In order to determine GR processes which occur in active 
galactic nuclei (AGN) an inverse problem of reconstructing the accreting system 
parameters from the line profile has to be solved quickly. In this paper we present 
a numerical approximation of the emission line Fe with analytical functions. 
The approximation is accomplished for a range of the disk radial coordinate r and 
the angle 9 between line of sight and perpendicular to the disk and allows one to 
decrease computing time by 10^ times in certain astrophysical problems taking into 
account all GR effects. The approximation results are available in the Internet at 
http : //www . iki . rssi . ru/people/repin/approx. 
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PACS: 98.54.Cm, 98.35. Mp 

1. INTRODUCTION 

In the past decade quite extensive X-ray observations of Seyfert galaxies have been carried 
out. In a considerable body of cases the wide emission line of Fe {Eq = 6.4 KeV) 
with a peculiar two-peak profile [l]-[20] is observed in the active galactic nuclei (AGN) of 
these galaxies. Maxima of the line are of different height, and a long red wing may stretch 
up to ~ 3 KeV. The Doppler line width corresponds to matter velocities of tens of 
thousand kilometers per second reaching v ~ 10^ km/s for the Seyfert galaxy MCG-6-30- 
15 [2] and v — 4.8 ■ 10^ km/s for MCG-5-23-16 [6]. Today the most reasonable explanation 
seems the one which indicates that the Fe line emerges in inner regions of an accretion 
disk (r ~ 1 ^ Avg) where GR effects dominate. The Fe Ka line is also observed in X-ray 
binaries and /xQSO's. A narrow Fe line may also be a part of a wide line, which cannot 
be distinguished from the background. 

In order to determine GR processes which occur in active galactic nuclei (AGN) an inverse 
problem of reconstructing the accreting system parameters from the line profile has to be 
solved quickly. To get the line profile one has to solve numerically GR equations of motion 
for photons with a wide range of initial conditions. These computations are very time- 
consuming. They have already been carried out before [21]- [35], but an amount of effort is so 
unusual that it is necessary to simplify essentially the computation procedure. Among other 
things, one has to decrease the computing time and the amount of computations. In this 
paper we present an approximation of the Fe line profile with analytical functions. The 
use of these functions offers a decrease of the computing time by 10"^ - 10® times depending on 
a line parameters. This procedure essentially simplifies the solution of the inverse problem 
of reconstructing AGN parameters from observational data. 

In Sec. 2 we consider the procedure of calculating the Fe line profile within the GR 
framework. In Sec. 3 we propose a method of approximating the Fe line profile with 
analytical functions. Sec. 4 deals with a numerical method (genetic algorithm) of search for 
optimal approximation parameter values. In Sec. 5 we present the approximation results and 
estimate its quality. In Sec. 6 we discuss limitations of applying the obtained approximation 
to the astrophysical problems. The summary is made in Sec. 7. 
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2. CALCULATING THE LINE PROFILE IN GR FRAMEWORK 



Space-time around a rotating black hole is described with the Kerr metrics: 



r'' + a'' + ^ sin^ o] sin^ 9 + ^ sin^ 9 d(t>dt, 
where (t, r, 6*, 0) are the Boyer-Lindquist coordinates. The standard notations are introduced: 



p2 = ^2 _^ Q,2 cogS A = r2-rgr + a^ rg = 2GM, 

where G is the gravitational constant, a and M are angular momentum and mass of the 
black hole, respectively. Hereafter we assume almost extreme angular momentum for the 
black hole a ^ 0.9981M and G = c = 1. 

Equations of motion of free particles in the Kerr metrics are obtained by separation of 
variables in the Hamilton- Jacobi equation: [36 38] 

i ik dS dS 2 

where m is mass of a particle. For photons m = 0. The Fe Ka line profile registered by 
a distant observer is obtained by solving equations of photon motion in the Kerr metrics. 
After we made the set of equations dimensionless we have: 

^ = -a (asin^^ - + (r^ + a' - ^a) , (1) 

dr , , 

n, (2) 



da 
dri 

da 
d9_ 

da 



2r'+{a^-e-v)r + {a-0^ + V, (3) 



01, (4) 

d9i J e 

—— — cos 9 — o 
da V sin"* 9 



cos^ ( -\- -a^sin^ ) , (5) 



da 



where a is an affine variable, r] and ^ are constants that define a trajectory of a particle. 
They are expressed through quantities conserving on the trajectory: the particle energy at 
infinity E, the projection of angular momentum on the z-axis {9 — 0) and the Carter [36] 
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separation constant^ Q, viz. rj — Q/M'^E^ and ^ = L^/ME. In the set of equations (l)-(6) 
the angular momentum a, the coordinates t and r are measured in the units of the black 
hole mass. Respectively, A measured in the units is related to the dimensionless r and 
a as follows: A = — 2r + a^. The extra variables ri and 9i are introduced for the set not 
to have singularities, and their physical meaning is not important to us. Two first integrals 
of the system 



are used to control the calculation accuracy and avoid accumulation of integrational errors. 
Namely, quantities ei and £2 have to be smaller than 10~^ in the end of the trajectory. The 
method of solving equations (l)-(6), results of simulations, as well as the derivation of the 
equations are available in papers [28, 39, 40]. 

To solve the set (l)-(6) numerically one has to set initial conditions. We assume that 
particles in the disk move along circular orbits. The disk is in an equatorial plane, optically 
thick and radiates monochromatic photons with the energy Eq isotropically in the laboratory 
reference frame comoving with the disk. We assume the disk to be opaque, i.e. the disk 
absorbs a photon crossing the disk plane. In other words, a photon emitted from one side 
of the disk cannot be registered from the other side. A profile of the line radiated by a thin 
ring with r — rQ = const is obtained by setting initial conditions at r = ro and collecting at 
infinity the photons, which come to a distant observer in the direction 9. A profile of the 
line radiated by the entire disk can be obtained by integrating the intensity with respect to 
the radial coordinate. 

Thus, given the disk and the direction 9, the line profile normalized at the unity in the 
maximum is the dependance of the intensity on the photon energy at infinity. The energy is 
measured in the units of the laboratory energy Eq: 




(a^ - - - 2 [{a -^f + ri\r + a^r] = 0, 



(7) 




(8) 



(9) 



where 



E 



(10) 




g = p2 _^ cos^ 9 



( 



— a^M^E^ ) , where is the ^-component of 4-momentum. 
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Figure 1. Typical Fe Ka line profile. Radial coordinate value r = 4r(,, angle between line of sight 
and perpendicular to the disk 9 = 30°. Photon energy in reference frame comoving with the disk is 
taken for unity. Small-scale fluctuations are due to statistical reasons. 




I I ' — i i i 1' I I 

0.5 0,6 0.7 0,8 0.9 1-0 1.1 1.2 

Energy 

Figure 2. Plot of function (12) witli ai = 0.011181, ai = 1.73588, 02 = 0.000194, 62 = 85.16571, 
a2 = 0.731736, Xi = 0.6711022, X2 = 1.005738. Parameter values correspond to approximation of 
the curve shown in Fig. 1. In this Fig. normalization is not taken into account. 

A typical Fe line profile registered by a distant observer with 9 — 30° from a ring of 
the radius Tq = 4rg is shown in Fig. 1. 

3. LINE APPROXIMATION 



If we take a look at the typical Fe line profile (see Fig. 1) we notice that each maximum 
has both sharp and gentle slope. Therefore, it is necessary to find an analytical function 
with this property. 
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As it is well-known from mathematical analysis, the function f{x) — exp{—a/x) tends 
to zero as X — > 0, moreover, it tends faster than any power-law function does. Besides, 
/(x) ^ 1 as x ^ oo. The power-law function g{x) = x~"-, a > 0, on the contrary, tends to 
infinity as a; ^ 0, and to zero as a; — > oo. Then their product 

y{x)^f{x)-g{x)^e-x.^[-^-^ (11) 

has all the required properties. On one hand, if a is large enough the latter function has a 
sharp slope (since the function behavior is exp(— a/x) if a is large), on the other hand the 
function has a gentle power-law slope (since exp(— a/x) ~ 1 if x is large). 

Writing the similar terms for each of the two maxima of the line we obtain an approxi- 
mation formula: 

, . X — Xi 1 X'2. — X ho , , 

y(-)-^ ■j^^^x^,+- -(Sj^- ■^'<^<^- (12) 

which comprises seven parameters. The parameters a\ and 02 determine sharpness of the line 
maxima. Xi and X2 determine approximate positions of both maxima, ai and 02 determine 
the function behavior between the maxima, and 62 sets the relative heights of the maxima. 
The coefficient hi may be set equal to unity. An example of a non-normalized plot of the 
function is shown in Fig. 2. 



4. NUMERICAL METHOD. GENETIC ALGORITHM. 

The approximation quality can be estimated by calculating a sum of squared deviations 
of the function (12) from a line profile, which is obtained by solving the set (l)-(6) and 
shown in Fig. 1. In other words, to estimate the approximation quality it is required to find 
a minimum of the function: 

(13) 

in the 7-dimensional parameter space (oi, Xi, ai, 02, X2, 62, ol2)- Here xi are selected points 
of the normalized energy (10), y{xi) is the approximation from (12), Umax is the maximal 
value from numbers y(xj), Ko,{xi) is the solution (9) obtained from the set (l)-(6), and 
n\ and 77-2 are the minimal and maximal index values when Ka{xi) 7^ (one can see in 
Fig. 1 that on the right and left sides the plot has zero segments; accurate meaning of the 
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variables rii and 712 is to be cleared below). In such cases a domain of the 7D space is usually 
covered with a grid and the function F values are calculated in nodes of the grid. However, 
it is practically impossible to accomplish the task with exhaustive search. For minimally 
acceptable approximation accuracy the grid should be fine enough: from to 2 with step 
0.01 in (ai, 02), from to 2.5 with step 0.01 in (ai, a2), from to 200 with step 1 in 62- The 
values (Xi,X2) are known approximately (they are positions of the maxima), however, even 
close to these values a grid with at least 100 nodes is required. Thus, for the acceptable 
approximation the whole number of the grid nodes (possible parameter combinations) is over 
10^^, which is beyond computing capability. 

For multidimensional problems like this one may apply the genetic algorithm, efficiency 
of which increases with an increase of number of dimensions. For the algorithm to work 
efficiently just continuity and relative smoothness is required. Using the algorithm does not 
guarantee that we find the exact minimum, however, this algorithm finds this minimum with 
quite high probability. 

The idea of the algorithm is taken from biology: the fittest survive, the weakest die. 
The domain of the multidimensional space is first covered with a sufficiently fine grid, and 
nodes of the grid are encoded with consequent binary natural numbers in each coordinate. 
This binary number with a certain number of zeros and units (bits) is called a gene. (e.g. 
00101101). Each coordinate gets its own number of bits. For instance, if the grid has 256 
nodes with respect to some dimension, the genes corresponding to this dimension contain 8 
bits (one of the genes, viz. the 46th, is mentioned above). The number of genes for each 
node is the same as the number of the space dimensions. A chromosome of a certain node 
is a system of the genes which are written consequently (in the parameters' order) without 
spaces (e.g. if 010 and 1001 are genes then 0101001 is a chromosome). The number of bits 
in a chromosome is equal to the sum of bits of the constituent genes. Value of a chromosome 
is defined as the value of the function —F in a node corresponding to this chromosome. 

The essence of the algorithm is the technique with which the chromosome with maximal 
value, i.e. the minimum of the function F, is searched. At the initial instant some number of 
chromosomes is randomly selected in the multidimensional space and put in descending order 
of their values. Then the operation of crossover is performed. To do this pairs of chromosomes 
are randomly selected; the higher a chromosome value the higher the probability to crossover. 
Crossover is an operation of cutting each of the two chromosomes at a random point and 
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exchanging the cut parts. For the newly born chromosomes their values are calculated. Then 
the entire population, parents and their off-spring, are put in descending order again, and 
the weakest chromosomes are omitted so that the number of chromosomes remains the same. 
After several generations the fittest chromosome corresponds with high probability to the 
minimum of the function F. 

Note that just like in biology chromosomes can mutate. Mutation is a random inversion 
of a bit in a chromosome. To fulfil the procedure we make the chromosome mutate after 
crossover, but before calculating its value. Namely, consequently scanning all its bits we 
invert each of them with very small probability (0.001 or smaller). If a generation contains 
200 chromosomes and each is 20 bits long then after mutation 4 bits (of 20 • 200 = 4000) will 
be inverted. This procedure allows one to refresh the generation and speed up the minimum 
search. See e.g. [44, 45] for detailed description of the genetic algorithm. 

The function (13) is calculated at n fixed points Xi uniformly distributed in logarithmic 
scale. We should first choose an interval between the points 

W ^-{IgEn - IgEL) ^- - Ig^, (14) 
n n El 

where El is the minimal value of the energy, Eh is the maximal value of the energy, n is 
number of the intervals, and fix points of the interval separation 

\gik^\gEL + kW, A; = 0,l,...,n. (15) 

The points Xj, in which we calculate values of the function F[x), are in the middle between 
the separation points: 

x,^, = ^^^^ = ^-^^ ■ 10^^^-+^^, = 0, 1, . . . , n - 1. (16) 

The values of the variables El and Ejj should be chosen in such a way that the spectral line 
were within the segment [El.Eh]- Even better if the segment covers the zero segments as 
well. To calculate values of the function F (13) not all values Xi are required, just the ones 
where Ka{xi) > 0. In other words, if rii and n2 stand for the minimal and maximal values 
of k where Ka{xi) > we use them in (13) as summation limits. 

Thus, the function values are calculated only at the aforementioned points Xi, and we let 
F{x) = as x < Xni and x > x^j. To perform the calculations the function y{x) should be 
normalized by dividing by the maximal value at the points Xi. The F(x) function behavior 



9 



between the points xi is not considered. Among other things, one cannot guarantee that 
the functions F{x) and Ka {x) have the maximum at the same point. Note that the interval 
width W in (14) also determines the precision of the function value localization of F{x) 
and Ka{x) on the x-axis. In other words, the values of these functions are considered to be 
constant between the points and ^k+i- 

The minimum of the function (13) was searched at several stages. At first, the minimum 
of F{x) was found 10 times with the constant parameter intervals on a quite rough grid 
(there were 1000 chromosomes and 80 generations). The best result of the ten was used to 
center the domain in the 7-dimensional space. After that the minimum search was performed 
again, but on a finer grid. Then the process was executed two more times. The result of 
the last iteration was accepted as the minimum of F{x). Note that such a cumbersome 
procedure is required, because the function F{x) has a number of accessory minima which 
are hard to distinguish from the main minimum (this is the reason why gradient methods are 
not applicable here). Using the genetic algorithm does not guarantee that we find the "very" 
main maximum. However, from our experience we can say that even if in some cases the 
genetics was able to find only an accessory maximum it yielded an adequate approximation. 

5. SIMULATION RESULTS 

In Fig. 3 the most typical results of the approximation for two sets of initial parameters 
r = 4.5 Tf,, 9 — 30° and r = 3.4 r^, 9 — 60° are shown. The initial theoretical curve 
obtained by solving numerically the set (l)-(6) is plotted with a thin line, and the result of 
its approximation is plotted with a thick line. Theoretical curves in this and other cases are 
wavy, but these fluctuations are statistical [40] (shot noise) and not physical. From Fig. 3 
one can see that the theoretical and approximation lines are quite close to each other and 
cannot be distinguished with a naked eye. 

The best approximation results are impossible to be shown on the plot at reasonable scale, 
because the discrepancy between the curves is about a thin line thick. One can imagine that 
one of these plots (with two curves) is shown in Fig. 1 or 2, but because of the tiny discrepancy 
between the curves we cannot see them separately. For r = 4.2 and 9 = 60° the curves 
deviate from each other at any point less than 0.5% from the maximal value. 

The results of approximating the line Fe K^, for the fixed value of the radial coordinate 
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Figure 3. Typical result of approximating the line Fe K^- The initial theoretical curve is plotted 
with a thin line, and the result of its approximation is plotted with a bold line. The values of the 
radial coordinate r and an angle 6 between line of sight and perpendicular to the disk are written on 

each plot. 
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Figure 4. Approximation results for the fixed value of the radial coordinate r = 3.4 and 
different values of an angle between Hne of sight and perpendicular to the disk from 9 = 10° to 
9 = 80°. X-axis scale is different for different plots. 

and different values of the disk inclination angle are shown in Fig. 4. One of the plots 
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Figure 5. Approximation results for the fixed value 6 = 40° of an angle between line of sight and 
perpendicular to the disk and different values of the radial coordinate from r = 0.8 rg to r = 10 r^. 

X-axis scale is different for different plots. 

with 9 = 60° shown in Fig. 3 can also be added herein. As one can see from the plots the 
approximation is adequate for all disk inclination angles up to 6' = 80°. In the latter case 
{9 — 80°) the proposed approximation model may become inadequate, because the lensing 
effects play a big role, and one or two extra maxima [46] may appear at some values of 
r. However, it is possible to apply the approximations for practical problems which do not 
require high accuracy. If one wants to obtain a more accurate approximation for curves with 
9 = 80° one can add to the expression (12) a quadratic function with a negative leading 
coefficient. This extra term will describe an extra detail in the spectrum (a detail similar to 
that in Fig. 4 for ^ = 80° close to E — 1.15). For 9 > 85° the proposed model is inadequate. 

In Fig. 5 the approximation results for the fixed value of the disk inclination angle 9 — 40° 
and different values of the radial coordinate r are shown. As one can see from the plots the 
approximation is adequate for a wide range of the radial coordinate values. Note that the 
approximation remains adequate at the very boundaries of the interval, i.e. for r = 0.8 Vg and 
r = 10 r^, and some statistical fluctuations at r = 0.8 r^, do not make it worse. Extrapolating 
the results we can guess that even for r > lOr^ the approximation (12) is quite reliable. 



12 




□ 1,4-1,6 

■ 1,2-1,4 

□ 1-1,2 

■ 0,S-1 

□ 0,6-0,8 

□ 0,4-0,6 

■ 0,2-0,4 

□ 0-0,2 



Figure 6. Approximation quality. Root-mean-square zi{r,9) in significant part of the curve 
measured in percentage, x- and y-axes are variables r and 0, respectively. 

In Tab. 1 values of the parameters oi, ai, a2, 62, «2, -^1, -^2 are given for some values 
of the radial coordinate r and the disk inclination angle 6. Some technical information, for 
which the parameter values were obtained, are given in Tab. 2. The detailed information on 
all the parameter values for 0.7 < r < lOr^ and 10° < 9 < 80° is available in the Internet 
at http : //www . iki . rssi.ru/people/repin/approx. 

To estimate the approximation quality beside the function (13) we also used 



zi {r,9) 



E 



Vmax 



(Xi 



(17) 



^2 - + 1 

where ni and n2 are the numbers of the first and last intervals in a significant part of the 
curve, respectively. This value is a root-mean-squared per interval in a significant part of 
the curve. The term "significant" means that we consider only the intervals with at least 
one photon. For instance, for r = 3.4 and 9 = 20° the calculations were carried out for 
El = 0.2, Eh = 1.4, n = 150. The line itself occupies the interval 0.66 < Xi < 0.92. The 
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statement made above means that while calculating the root-mean-square the intervals with 
Xi < 0.66 and Xi > 0.92 were omitted. 

It is convenient to measure the function Zi (r, 9) in percentage. The plot of this function 
is shown in Fig. 6. The values 0.2 - 0.5% are typical. The best approximation results are 
for l.bvg < r < Gvg. For higher and lower values of r the approximation is some worse, but 
even though the deviations higher than 0.7% are quite rare. For stand-alone points with 
r < Tg and r > 6.5 Vg a more accurate approximation will be found some day. 

To estimate the approximation quality one can use an extra criterion which describes the 
maximal discrepancy between theoretical and approximation curves 



Z2 (r, 9) — max 



i^l,...,n. (18) 



Umax 

and measured in percentage of the maximal value of a curve. The plot of the value is 
shown in Fig. 7. As one can see from the plot, the maximal discrepancy is mainly 6-8%, 
and discrepancies higher than 10% are quite rare. Besides, this criterion has no explicit 
dependance on variables r and 9. Intervals with maximal discrepancies are most often 
located close to the blue (higher) maximum from its left side where the curve drops sharply. 
Usually it is one or two points. Thus, this criterion characterizes "the worst" point of the 
curve. 



6. DISCUSSION 

As it has been noticed above, the function F{x) values can be taken only at the points 
Xi defined by the relations (14) and (16), because between these points the function F{x) 
values are not defined. Since it also applies to the points Xi close to the maxima the used 
approach does not guarantee that the function F{x) reaches its maximum at the same point 
as Ka (x). When using the approximation (12) (e.g. to model a profile of the line radiated 
by the entire disk [41-43]) one should first calculate intensities at the points Xi for several 
thin rings and then sum up the results with a decrease of number of intervals in the resulting 
curve (hence, make it rougher). This procedure allows one to obtain a quite smooth curve 
which can be used to interpret observational data. 

In Fig. 8 examples of the approximations with the highest deviations from the theoretical 
curves are shown. Although these cases are rare, we find it important to explain reasons 
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Figure 7. Approximation quality. Plot of the maximal discrepancy (in percentage) between curves 
is shown, x- and y-axes are variables r and 6, respectively. 

which caused them. One of the reasons is improper photon energy distribution in higher 
(blue) and lower (red) maxima. As a result, a maximum consists of two values with almost 
the same height which is not suitable for the approximation (12). This situation occurs for 
the second (red maximum) and third (blue maximum) plots in Fig. 8. This problem can be 
solved by changing the location of the points Xj close to the maxima which, in turn, can be 
achieved by changing values of El, Eh and n. In this way one can obtain approximations 
of higher quality for the curves shown in Fig. 8, but this is quite time-consuming. The 
other reason appears at the small values r < 1.2 r^. Here in a quite wide range of the disk 
inclination angle ^, as a matter of fact, a curve has just one distinct maximum and the 
approximation (12) does not work. 
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Figure 8. Examples of approximations with the highest deviation from theoretical curves. Values 
of radial coordinate r and angle 6 between line of sight and perpendicular to the disk are written on 

each plot. 

7. CONCLUSIONS 

We obtained an analytical approximation of the iron emission line profile in AGN, which 
allows one to determine physical parameters of the matter of the accreting disk around a 
black hole by the line profile. 

According to the criterion (17) in the intervals 1.3 < r < 2.5 with respect to the 
radial coordinate and 10° < 9 < 80° with respect to the angle 9 as well as 3.5 < r < Gr^, 
20° < 9 < 80° the approximation (12) yields the accuracy of 0.8% excluding several points 
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with higher deviation, but still lower than 1.4%. In the interval 0.7 < r < 1.3 for all 
values of 9 one has the accuracy of 1.2% excluding two points with 1.6%. 

According to the criterion of the maximal discrepancy (18) the approximation (12) yields 
the accuracy from 2%. to 12%. in the entire investigated interval 0.7 < r < 10 r^. No 
explicit dependance of the approximation accuracy on r and 6 was noticed. As an exception 
there are less than ten points (mainly at r < Avg) in which the discrepancy amounts to 18%. 

The application of the approximation (12) in practical astrophysical problems results in 
a decrease of computing time by 10^ — 10^ times. 

One of us (SVR) expresses his gratitude to Prof. E.V. Starostenko, Dr. R.E. Beresneva and 
Dr. O.X. Sumenkova for the possibility to work intensively and fruitfully on the considered 
problem and to Prof. A.F. Zakharov for fruitful discussions. VNL and VNS are grateful 
to Russian Fond of Basic Research for partial support of this paper, Grants 04-02-17444, 
07-02-00886. VNS thanks UNK FIAN for support. 
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r 9 ai ai 02 62 a2 Xi X2 



1.0 80 0.046261 0.15376 0.024139 97.82825 1.304885 0.1474317 1.416187 
3.4 10 0.011971 1.7246 0.000971 37.33291 0.766634 0.7158835 0.8389752 
3.4 20 0.010876 1.72776 0.000365 84.92977 0.672685 0.6712679 0.9066637 
3.4 30 0.008042 1.41 0.00064 35.70628 0.701063 0.6354619 0.9801457 
3.4 40 0.011937 1.569 0.000698 62.20945 0.730865 0.6015416 1.059468 
3.4 50 0.011593 1.7134 0.001139 37.30477 0.768584 0.7159498 0.8390415 
3.4 60 0.006658 1.4788 0.001595 82.66419 0.845563 0.5537415 1.222009 
3.4 70 0.003384 1.4204 0.005798 84.84075 0.939793 0.5395592 1.287379 

3.4 80 0.000955 1.15504 0.004171 43.38763 0.90118 0.532406 1.321832 

3.5 10 0.008442 1.67216 0.000816 56.68913 0.731884 0.7252888 0.8497934 
3.5 20 0.005904 1.40092 0.001682 29.44844 0.73757 0.680184 0.9187329 
3.5 30 0.012463 1.65552 0.002242 65.77657 0.756707 0.6386765 0.9931731 
3.5 40 0.014688 1.66824 0.011318 51.58594 0.959569 0.6063802 1.073101 
3.5 50 0.002631 1.28648 0.00097 40.9329 0.802594 0.5832343 1.145212 
3.5 60 0.007254 1.44316 0.002104 63.847 0.853285 0.5609719 1.222076 
3.5 70 0.006386 1.5494 0.013498 97.622 1.032985 0.5466044 1.286815 
3.5 80 0.001006 1.33604 0.00075 77.44388 0.890169 0.5394267 1.320739 

4.1 10 0.00591 1.517 0.000592 49.33127 0.643282 0.7605556 0.8719551 
4.1 20 0.004271 1.29656 0.00042 31.59231 0.61765 0.7180462 0.9353039 
4.1 30 0.006604 1.37288 0.000688 36.44845 0.648607 0.6786068 1.00382 
4.1 40 0.004174 1.37864 0.001503 62.75788 0.692406 0.6498227 1.075424 
4.1 50 0.001696 1.17712 0.000427 39.10009 0.685844 0.6245808 1.1435 
4.1 60 0.012402 1.63436 0.002646 89.25638 0.771792 0.5985646 1.20719 
6.5 60 0.019516 1.66296 0.0003 34.00781 0.859408 0.6902458 1.191818 



Table 1. Approximation results for some values of radial coordinate r and angle 6 between line of 

sight and perpendicular to the disk. 
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Table 2. Technical information for the same values of radial coordinate r and angle 6 between line 

of sight and perpendicular to the disk as in Tab.l. 
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